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Abstract 

We extend to finite temperature the time-dependent effective theory for the Gold- 
stone field (the phase of the pair field) 9 which is appropriate for a superfluid containing 
one species of fermions with s-wave interactions, described by the BCS Lagrangian. We 
show that, when Landau damping is neglected, the effective theory can be written as a 
local time-dependent non-linear Schrodinger Lagrangian (TDNLSL) which preserves the 
Galilean invariance of the zero temperature effective theory and is identified with the 
superfluid component. We then calculate the relevant Landau terms which are non-local 
and which destroy the Galilean invariance. We show that the retarded ^-propagator (in 
momentum space) can be well represented by two poles in the lower-half frequency plane, 
describing damping with a predicted temperature, frequency and momentum dependence. 
It is argued that the real parts of the Landau terms can be approximately interpreted as 
contributing to the normal fluid component. 
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1 Introduction 



The form of the effective (Ginzburg-Landau) Lagrangian for time- dependent superconducting 
phenomena continues to attract attention. At T = Schakel using derivative expansion 
techniques 0, obtained a time-dependent effective Lagrangian for the Goldstone field ( which is 
the phase 9 of the pair field A ), starting from a simple BCS Lagrangian. Interestingly, Schakel's 
effective Lagrangian is the same as one that had been proposed earlier || by assuming that the 
momentum density in a Galilean invariant system is proportional to a conserved current. In [|J 
(see also ||) the Lagrangian of Jj]] and |3j was shown to be equivalent, up to a certain order in 
derivatives, to a time-dependent non- linear Schrodinger Lagrangian (TDNLSL); in particular, 
the role of Galilean invariance was emphasized. More recently, de Palo et at. || have given 
a more general derivation of the TDNLSL at T = 0, by including in the effective action from 
the outset fluctuations in the density field p (conjugate to 9), as well as those in A. This 
formulation allows a unified description of the weak-coupling (BCS) and the strong-coupling 
(Bose) limits as well as of the cross-over between them |J. 

In all the above cases, at T = 0, electromagnetic interactions can be straightforwardly 
included by the "minimal coupling" prescription, since the effective Lagrangian is local. As is 
well-known [[/]], the q 2 = pole of the #-mode is then pushed to finite frequencies, becoming 
the plasma mode. 

At finite temperature, in the non-static case, the situation is complicated by the presence 
of Landau damping terms, which appear to prevent any local time-dependent Lagrangian for 
9, since they have singularities at go = ±i>iH q| ( see Section 6 below), and the resulting cut 
structure implies that there is no Taylor series expansion about the origin in momentum space. 
This problem was first pointed out by Abrahams and Tsuneto ||. It poses a fundamental 
difficulty for the inclusion of electromagnetism, since minimal coupling cannot be used unless 
the effective Lagrangian is local. Note that the singularities, and hence this difficulty, disappear 
in the static case, so that there is, of course, nothing wrong with the familiar static Ginzburg- 
Landau theory at T ^ 0. 

Adopting a modern field-theoretical approach to the problem of deriving an effective La- 
grangian, Stoof H neglected the Landau damping terms entirely and constructed a time- 
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dependent effective theory by writing the pair field as A(x, t) = \Aq\ + A'(x, t), where |A, 
is the value of A which minimizes the effective action, and A'/|Ao| was assumed to be varying 



slowly in space and time. In [I0| an attempt was made to extend Stoof's approach so as to 
include Landau damping, retaining only terms quadratic in A'. It was argued that the trou- 
blesome (Landau) cut structure could be replaced - in the ^-propagator - by two simple poles 
at complex energy. This led to an approximate local effective action for 9 of the form 

J d*x [0 2 - [a(T) - ib(T)]v 2 F (V9) 2 ] (1) 

where a(T), b(T) are real functions. If ([!]) were correct, minimal coupling could, as usual, be 
employed to include electromagnetism. 

Unfortunately ([!]) is not correct. The imaginary part of the Goldstone mode propagator, Sq, 
is an odd function of the frequency and will not contribute to the total effective action as in (jl]) 
(it does, however, enter the equation of motion for 9). The correct approximate two-pole form 
of Sq will be given in Section 7 below, and it will be such that Sq 1 is not local in coordinate 
space. Thus we are back with the problem of how to include electromagnetism in the non-static 
case for T ^ 0. 

Since we do not at present know of any answer to this problem, all our subsequent analysis 
will generally refer to BCS-type pairing of neutral fermions, i.e. to an uncharged fermionic 
superfluid. Helium-3 is such a fluid but it condenses into a p-wave state, not the s-wave BCS 
state studied here. It would be straightforward, and interesting, to generalize the present 
approach to p-wave pairing (a beginning was made by Van Weert (ll]]). However, it is likely 



that trapped condensates of neutral fermionic fluids will be produced before very long (see for 
example [0). Our analysis would be relevant to such a case, when re-done for a finite system. 

Apart from the error in Sq, the approach leading to (P is unsatisfactory for another reason, 
namely it fails to reduce, at T = 0, to the correct Galilean- invariant result of [1], ^|,|J. The 
problem can be traced to the step A = |A | + A': as was pointed out in M, this ansatz 
destroys Galilean invariance and fails to allow for large fluctuations in 9, two weaknesses which 



are in fact related. In [10], the further approximations A fa |A| (1 + %9) and |A| = constant 
were made, which led to the form (|IJ) after retaining only quadratic terms in Im(A). What is 
required, rather, is the finite temperature version of 0] (or even better of ||), where we write 
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A = |A| e te , ignore fluctuations in |A| and consider a systematic expansion in the derivatives 
of 9, without having to assume that 9 itself is small. 

Our aim in this paper is to provide the finite temperature extension of |4[] and arrive at 
a temperature and time-dependent effective Lagrangian for 9, L{9). In Section 2 we set up 
the effective action formalism and derive the usual gap equation. In Section 3 we evaluate the 
first term in the expansion of the effective action in powers of derivatives and in Section 4 we 
calculate the next term which includes the Landau damping contribution. In Section 5 we show 
that, if the Landau terms are neglected, Galilean invariance is maintained (at least up to the 
order in derivatives to which we are working) and L(9) has exactly the same form as in 
but with temperature dependent coefficients. This, in turn, implies that an effective non-linear 
Schrodinger theory can be obtained with a temperature-dependent potential term. 

Next, we turn our attention to the Landau terms, which are calculated in momentum space 
in Section 6. As expected, their form is similar to the expressions found in [JToJ] . Their presence 
prevents any representation of the effective action in terms of a local Lagrangian. In Section 7 
we consider how to obtain a useful approximate form for the equation of motion for 9, including 
the Landau terms. Noting (see above) that Ihi(S'g') is an odd function of frequency, we show 
that Sq can be well represented by two poles in the lower half plane of the complex frequency 
variable. This corrects the form proposed in and invalidates (0). We then consider the 
real parts of the Landau terms, neglecting the imaginary parts. We propose that they can be 
approximated by their value at the (undamped) poles of Sq - an "on-shell" approximation. 
The complete effective action then becomes local, though the Landau terms break Galilean 
invariance. We propose a physical interpretation of the resulting equation of motion, in which 
the real parts of the Landau terms represent a contribution to the normal fluid component. 
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2 The Effective Action and the Gap Equation 

We begin with a brief review of the formalism. The BCS Lagrangian, for s-wave pairing and 
in the absence of external fields, is 



L = COX*^ + — + fJL)^ ff {x) + g^{x)^l{x)'ij} i {x)^{x) 
CT=1 Zm 



(2) 



where ij) a (x) describes electrons with spin a = (f, |), // = ^ is the Fermi energy, x = (t, x) 
and g > 0. We define the Nambu spinor 



and introduce the auxiliary "pair" fields A (re) and A* (a;) so as to write the partition function 

as 



Z = J Dip^Dip J DA*DA exp % J d 4 x (^G~^ 



(3) 



where 



G" 1 

— ' T.:r 



\ 



We now do a gauge transformation on the electron field writing ip a = exp f^y^) X<j( x )) 
where 9 is the phase of A and \ a are new anticommuting fields. The integrals over \ an d X 
are then performed, leading to 



Z = J DA*DA exp(iS eff [A,A j ]) 
where the effective action S e g is given by 



(4) 



S eff [A, A*} = -zTr log (CT 1 - E) - / d 4 x 



IAI 



— zlr iogG + zlr 2^ I a x- 



n=l 



n 



(5) 
(6) 



with 
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./ 




\ 2 ' 



1 2 8m 1 ' 



Note that the effect of the gauge transformation has been to separate the dependence on |A| 
which is in G~ l from that on 9 which is in S. We point out that Det and Tr refer to both 
the functional space and the internal (Nambu) space. The operators po, P and x are Fock- 
Schwinger Proper Time operators. 

The ground state is determined by taking |A| to be space-time independent and solving the 
mean field equation 

ss eff _ n ^ ■ f d4 P 1 _ i m 
~6\A\~ 9 J (27r)*pg-^- 1 (?) 

where E 2 = e 2 + |A| 2 and e = |— — ^ . The modulus |A| of the pair field is therefore the energy 
gap. (0) can be written in the usual form of the BCS gap equation at finite temperature ]7|, 13 



f d 3 p 1 BE 

q I - — r^r — ~ tanh = 1 (8) 

y J (2tt) 3 2E 2 V ; 

where we assume that the potential (the coupling constant) is attractive in a shell of width 2lu c 
centered at the Fermi surface and vanishes outside this shell. We note that in the weak-coupling 
limit uj c ^> |A(0)|, where |A(0)| refers to the value of the gap within the shell. In this limit we 
can consider |A| to be a time and space independent constant whose value is 

|A(T)| ~ |A(0)| ~ 1.75k B T c (9) 

for 



b = T/T c G [0.1, 0.6]. (10) 

Thus the effective theory only involves the single degree of freedom 9. If we expand each 
term of the infinite series in (||]) in powers of momentum and frequency and restrict ourselves to 



small values of these quantities, it has been proved in ]T],|J that at T = the trace of the first 
term contributes a term proportional to [9 + (V9) 2 /4m] to L{9), and the trace of the second 
term contributes a term proportional to [9 + (V9) 2 /Am\ 2 . As emphasised in Q, the combination 
[9 + (V#) 2 /4m] is Galilean invariant and it is this feature which allows L{9) to be rewritten as 
a TDNLSL at zero temperature. We shall now extend this calculation to finite temperature 
and explore to what extent an effective Lagrangian is still derivable and - if so - whether it can 
be written TDNLSL. 



3 The First Trace 

In this section we calculate the first (n = 1) term of the sum appearing in (|6]) : 

Tr[GS] = / P )n] / d 4 x \[9 + &$L\ (11) 

where a term odd in p has vanished in the last step, since we integrate over p from — oo to 
+oo. By tr we mean the trace only in the internal (Nambu) space. We apply the Matsubara 
formalism to evaluate the first integral at finite temperature. We Wick-rotate via po — > ipoE 
and then replace the integral over p QE by a sum over discrete fermion frequencies via 



oo 



oo 



^T-WPob) "4 E f^n) (12) 

27F P n=-oo 



where z n = (2n + l)ir/(3, and n is integer. Next, following |14], we write 



^J2f( iz n) = ^J2fz( z n) = -T- I dz f E (z) tan (^) (13) 

where the contour Cr encloses the real axis in a counter-clockwise sense. (JT3|) is valid provided 
that /e(z) has no singularities on the real axis, which is our case. It is convenient to rewrite 
the above equation in such a way as to separate out the (known) T = contribution : 

— / dz f E (z) tan(— ) =L1 + L2 + L3 (14) 
An Jc R 2 

with 
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LI 
L2 
L3 



dz f E (z) n(-iz) 
dz f E {z) n(iz) 



-oo 

2lX J — oo+ir/ 
2 poo—irj 

2lX J —oo—iri 



(15) 



where r\ is an infinitesimal positive quantity and n(x) 



i 



l+exp fix 



is the Fermi-Dirac distribution 



which has been introduced in our calculation through the use of the identities tan (4r^ 
-z[l - 2n(iz)\ = i[l - 2n{-iz)\. 

Applying the above formalism, we find that the total contribution of the first term of the 
series for S e ft is 



iTr[G±] 



d 3 



P 



(27T) 



E 



[1 - 2n(E)} 



2 [ + Am 1 



(16) 



which has the same Galilean-invariant form as at T = but with a temperature dependent 
coefficient. In fact, the temperature dependence of the first trace is not important and can be 
neglected. The maximum value of the function n(E) for the range of temperatures of (|10"D and 
with the gap having the value given in (RI) is n max ~ 0.05. Consequently 1 — 2n(E) ~ 1 and 



the coefficient of the Galilean invariant term [8 

d 3 p 



(V6>) 2 
4m 



IS 



M(T) 



(27T) 



E 



[1 - 2n{E)\ 



Po = M(0), 



(17) 



where po — VfI (3 71 " 2 ) is the electron density at T = 0. 



4 The Second Trace 

The trace of the second term (n = 2) of the infinite series in (^ can be written as: 

Tt[±G±G] = tr / ^ ^ E(g, p - |)G( Po - f , P - \)%(-q, P + ^)G(p + f , P + \) (18) 

where S(p, q ) = / c? 4 x exp (ip.x) E(x, q ). After performing the Wick rotation we have 
T r^AAAi f d 4 q E d 4 p E rr . f d 4 q E d 4 k E d 4 p E f d 4 q E d A k E d A r E d A p E 

TT[^G\ E --J ^ - 2zy (27r) 4 (27r) 4 (27r) 4 T 2 + _/ (27r) 4 (27r) 4 ( 27r )4 (2^ T » ^ 
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where 



Ti = 9 E {qE) E {-qE) A E ~ A E -\ 



BE[ (^ + £k ]+CE[ _(^)! + £k ] + ^ 0E P^ 



T 2 = 9 E (q E ) 9 E (k E ) 9 E (-q E - k E ) A E „ A EA 



(k q 



5m 



[B E + C E ] q 0E + D E ^- 

m 



T 3 = 9 E (k E )9 E (r E ) 9 E (q E - k E ) 9 E (-q E - r E ) A E „ A E _ 



-k 2 + k .q 

32m 2 



> 2 + r .q)[B E + C £ 



with 



2 i — ' Ae+ = A E (p E + A E - = A E (p E - — 



A e{pe) ■ , , 

POE + MP ) 

B E (poE,qoE,P , q ) = |A| 2 



C E (poE,qoE,p ,q 



-e(P - y)e(P + — ) + [Pqe ~ (— ) ) 



q_ 

2 



D e (poe, qoE, P , q ) = e(p + y)(p i? - ^) + e(p - y)bo£ + 
We note that the integrals over p E ah have the form 



9e 



dp t 



'OE 



A E _ A E+ g E (p E,qoE,p ,q 



9 = B,C,D. 



(20) 



'-oo 2tt 

We now proceed to the second step of the finite temperature formalism which is to discretise 
the momentum variables as follows 



P«e^pIe = ^^ , W2^oV2- (2m7 ° 



, j i (2/tt) s (2 S7 r) 

fas -> ^ = > r o£ -> ^ B = ^— ^ (21) 



and then calculate the sums over n which are contained in ([20]) by replacing them with integrals 
over z according to (p~3|) . The poles of tan(4^) are on the real axis and the sum of their residues 
is equal to the finite temperature sum which is to be calculated. We note that only A,B,C,D 
depend on the variable z and B,C,D have no poles. This means that we need to evaluate only 
the following integrals 

~g E = jf dz A E (z + M p + £.) Ae (z - M p - SL) g E (z, q^ p , q ) tan (^) (22) 
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where the function qe has no poles (analytical everywhere in the z-plane) and the A-functions 
have the poles 



2 ' 2 

where we have introduced the usual notation 



■^+iE + , -^-iE + 



E 4 



E(p + y) 



E. 



E(p - 



q 



We follow the same procedure as for the first trace. We note from (21) that exp (if3^jf- 



and then, following |f7|, pjp^], we Wick rotate the above expression back to real continuous 
variables go via q™E ~~ * ~iqo- The result can be written as g — §r + §l, where the regular part 
cjn has a well-defined expansion about the origin (go, q ) = (0, 0), but the Landau part g~i does 
not, as was first pointed out in j8|. For g = B,C,D we have 



with 



Br 
B L 
C R 

Dr 



IAI 



[1 - n(E + ) - n{Ej)\ 



1 



AE + E_ 

|A| 2 n{E + )-n{EJ) 
2E + E_ E + -E_+q 

Fi[l-n(£+)-n(£L)] 

2F n(E_)-n(E + ) 
2 E + -E_+q 

G x [1 - n(E+) - n(£L) 

n(E + ) -n(E_) 



+ 



1 



2G 



E + -E_+ q 



E + + E_-q E + +E_ + q 0j 



+ 



E.L.E- 



1 

4 



1 



E + + E_ - q E + + E_+ q 



E + + E_- q E + + E_+ q 



1 + 



E.L.E- 



G 1 



2E 4 



E, E_ 



After the analytic continuation, wherever we write g i n the denominators of B, C, D, we imply 
that there is an infinitesimal imaginary part, i 5, 5 — > + added to q , corresponding to retarded 
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boundary conditions |7|,|5|,|I5]]. We complete the Wick rotation back to continuous variables via 
k™ E — > — iko, r™ E — > — iro, and the second trace finally takes the form 



Tr[EGSG] 



+2z 



(2tt 



(27T) 



IP q 



2m? 



+ D (-i)go 



2m? 

p q 

m 



d4q d ' k e( q) e(k)0(-<i-k)^ k - cl+k2) 



(2vr) 4 (2tt) 4 
rf 3 p 



5m 



(2tt) 3 
d 4 g <i 4 fc <i 4 r ~ 



[5 + C](-<j) go+J D- P ' q 



m 



(2tt) 4 (2vr) 4 (2tt) 4 
(-k 2 + k .q )(r 2 + r .q 



0(fc) 0(r) %- fc) 6(-q-r 
d 3 p 



32m 2 



(27T) 



5 + C 



(23) 



5 The effective Lagrangian at T ^ 0, neglecting the Lan- 
dau terms 



The contribution of the first Trace to the effective action is given in (llq). When the Landau 
terms are neglected, only the functions Br, Cr and Dr appear in the expression (|23|) for the 
second Trace and they can all be Taylor-expanded about q = in powers of qo/ | A | and 
vf I q | / | A |, so long as these quantities are small. For Br and Cr the expansion takes the 
form 



B R {p,q) 
C R (p, q) 



A 



4£ 3 
| A | 2 

AE 3 



;i-2n(E)) + 0(q 2 ) (q.p) 2 ,g 2 ) 
;i-2n(£)) + 0(q 2 ,(q.p) 2 ,g 2 ). 



Dn(q), on the other hand, vanishes as q = 0, and we have 



(24) 
(25) 



D R (p,q) = O(g (q-p)). 



(26) 



When (j24|) - (|26 ) are inserted in (^) and the appropriate Fourier transforms performed, we 
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obtain a local effective Lagrangian for 6, involving powers and products of the derivatives 9 
and V0. 

We must now consider how this derivative expansion is to be ordered. Due to the preferred 
frame supplied by the heat bath, we do not expect Galilean invariance to be preserved in the 
effective Lagrangian at T 7^ 0. Nevertheless, we expect Galilean invariance to remain approxi- 
mately valid for T < T c . Accordingly we shall regard terms of the form and "(V#) 2 /m" as 
being of the same order - or equivalently, terms of the form "go#" and "q 2 # 2 " as the same order. 
It then follows that the leading order terms in the expansion of ( p3|) are given by retaining only 
the zeroth order contributions to Br, Cr and Dr from (23))-( p6|) , thanks to the explicit factors 



involving q and q already appearing in fl2"3"|). This is an important point. 



In evaluating these leading order terms, it is convenient to define 



(2tt) 3 2£ 3 
and to make use of the approximation |J 



d 3 r> I A I 2 

N ( T ) = I TTS-. hJr^-ME)), (27) 



2 

where 



N(T) « \n(0)R(T) (28) 



mpF 3/9, 

is the density of states at the Fermi energy and 



"<°> = # = ^ (29) 



nm= Idy 1 -^ (so) 



where y = e/ | A |. The leading contribution to the first term in (^) is then 

- i\N{T) J ^ ~e{q)~e{-q)ql = ~N(T) J d\ 6\x). (31) 

Evaluating the other terms in the same way, and including the expression ([TBD for the first 
Trace, we obtain the effective Lagrangian (neglecting Landau terms) 
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C R (9) = —M(T) 



6 + 



4m 



+ 



N(T) 



9 + 



(vey 

4m 



(32) 



Remarkably, Galilean invariance is still preserved by (^). The only difference between 
and C e jj{6) of |j],f| is the appearance of the T-dependent coefficients M(T) and N(T). In fact, 
as T -»• we note that M(T) -> p /2, while JV(T) -> iV(0), so the T = result of 0,| is easily 
regained. 

The Galilean invariance of (|32|) guarantees that it can be rewritten as a time-dependent 
non-linear Schrodinger Langrangian (TDNLSL), as explained in [|J. The results of can be 
immediately taken over to T ^ (always with the neglect of the Landau terms) by making the 
replacements po/2 — > M(T),N(0) — > N(T). The result is merely that the "potential" term in 
the equivalent TDNLSL becomes 



_ [p~2M(T)] 2 
K " 2iV(T) • 

Thus, always with neglect of Landau terms, we have shown that the effective theory for our 
neutral fermionic fluid is a TDNLSL, just as it is for a bosonic fluid. The inclusion of electro- 
magnetism by means of the "minimal coupling" prescription is straightforward. 

We now turn to the Landau terms which, as in | I0[] , are non-local and cannot be Taylor 
expanded about q = 0. 



6 The Landau terms 



The Landau terms B^, Cl and cannot be treated in the same way as the regular terms, 
because of their non-analyticity around the origin. However, since we are interested in the low 
frequency, long wavelength regime, we can follow [p,[l~4l and write 

e p.q 



E,-E_ 



(34) 



E m 

We now expand the numerators of these quantities in powers of momentum q and keep the 
leading order terms which are of first order in q. We find 



B, 



t; h ?o + w 

E m 



- 1 |A| 2 dn e 



p.q 



2E 2 dE E m 
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c L 

Dl 



L l 23 

E m 



1 -l 



+ q + iS 



2e 2 + |A| 2 dn e p.q 



~B Vqo + io 

hi m 



_1 dn ie 2 



2E 2 dE E m 
p.q 



(35) 



dE E 2 m 

If we look at (|23f) and recall that g = c/r + <Jl for g = B,C,D we realize that the Landau terms 
appear in the second Trace in the following integrals 



h 
h 
h 



d 3 p 

(2vr) 3 

d 3 p 

(2tt) 3 

d 3 p 

(2^ 



B, 



2m 2 



p q 



2m 2 



-*j9o- 



p q 



[B L + C L ](-t)q + D L ? q 



Substituting into Ji the expressions (]35f) for B L ,C L , D L we obtain 



n? 



(36) 



|p| 2 <i|p| rfn 
4tt 2 



where 



|p| 2 |q| 2 |A 

3m 2 

c/,s 



£ 2 



1 + f I.) % 2 



q Em 



(37) 



'-i s — a — z<5 sgn(e) e |p||q| 
The integral I s has a pole at a, if —1 < a < 1. The sign of e appears in the denominator of I s , 
because the infinitesimal imaginary part iS is multiplied by a whose sign depends only on the 
sign of e, the rest of the parameters in a being always positive. Performing the s-integration 
yields 



In I- -I + m sgn(e) 0(1 — \a\). 

1 + a 



(38) 



Similar expressions can be found for I 2 and I 3 . Applying now approximation (ESI) we find 



h 
h 

h 



^) (^JH 2 -{mc)+iHl 1 \c))ql/vl 
■ i 3p 



'2mvp 



{H?\c)+iH?\c)) 



(3), 



c)) 



(39) 
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where 



J 



i a i f , dn 



-\A 

\A\tic 

vf |q|' 

r(«) 



_ c y/y' 2 + 1 y + cVy 2 + 1 
2 



1 1/1 



y - cVy 2 + 1 



(40) 



We observe that, for |c| > 1, = , so that the cut in the functions is over the region 
|c| < 1. It is noteworthy that the function which multiplies the term in g 2 , is, apart from 
a factor of 4 and slight changes in notation, exactly the same as the corresponding function H 
of [Kj. 

Substituting the expressions which we have found for Ii,I 2 ,l3 back into the second trace, 
making approximation and Fourier transforming the 9 functions wherever this is possible, 
we can bring the effective action to the following form 



S P 



d x 



3p c 



5WI 



R 

2^ 



6 + 



(Vg 
4m 



-/ ^A<l)h-<l)^{HV{c) + iHl l \c)) 



(2tt) u f 

_a ^ s - to) «(*) S -(- 9 - , )290 M±i!) : L (J/ « (c) + imc)) 

d A q d 4 k d 4 r ~ 



(2tt) 4 (2tt) 4 (2tt) 4 
-k 2 + k.q) (r 2 + r.q) 1 



0(k)6(r)0(q- k)0(-q-r 
{H?\c)+iH?\c)) 



, ■ (41) 

Am Am vp J 

Consider now the effect of the Landau terms J and in (f4T|). The simplest of these is J, 
which is a (temperature-dependent) constant. Thus its contribution is still local but its presence 
does break Galilean invariance, since it disrupts the invariant combination U 8+(V8) 2 /Am" in the 
first term in square brackets in (|4lD . Actually, as far as this single term is concerned, Galilean 
invariance can be formally restored by modifying the first 9 term so as to read [1 + J\9 ; this 
does not affect the equations of motion. 
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The effect of the functions is - depending on their numerical magnitudes - potentially 
much more drastic. Most importantly, these functions cannot be expanded about q = 0, since 
they depend explicitly on the ratio c = qo/vp\q\. Thus the terms involving the if' n )'s in (|4T| ) 
cannot be represented as any Lagrangian local in coordinate space. A fortiori then, these terms 
also break Galilean invariance. They also contain imaginary parts, of course, corresponding to 
the physical " Cerenkov " radiation process in the heat bath P,|lT|. In the next section we 
discuss how the equation of motion for 9, which follows from fl4lf) , may be approximated by a 
more manageable form, which may be given a physical interpretation. 



7 Approximate equation of motion for 9 and its inter- 
pretation 



7.1 The ^-propagator 

Let us begin by considering the terms in (f4T|) which are quadratic in 9, from which we can read 
off the momentum-space propagator for 9. The first point to note is that, since is an odd 
function of q , it contributes nothing to the action S e g. To see the effect of damping we need 
to consider the equation of motion for 9. This has the form 



9 = S G (q , |q|) 



? 2 H${c), 



9 3 H?)(c) 



(42) 



where we have written the non-linear terms on the right hand side in symbolic form and where 
the ^-propagator is given by 



Sc(qo, |q|) = \ q 2 



At T = we have R = 2 and J 



c 2 R 1 + J 



2 3 



+ c 2 (HV>(c)+iHP(c 



(43) 



H\ = 0, so that Sg has poles at c = 



as 



expected for the Bogoliubov- Anderson mode fli~6|,|T7|j. 



Consider now what happens as T moves away from zero. ^ and J deviate slowly from 1 and 
respectively (see fig.|l]a, fig.|l|b ) while the quantities c 2 H^ and c 2 H^ are very small compared 
to the main term — | (see fig.^ja, fig.^b ). This suggests that it should be a good approximation 
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to replace the complicated c-dependence of the denominator in (|^) by a quadratic function, 
which amounts to replacing the cut structure of (|43[) by two poles - the same two that were 
present at T = but now at slightly different positions. In particular, the poles will now be 
shifted off the real axis, due to the presence of the function. It is important to note that 
Hjp is an even function of real c while iff is an odd function of real c. It is then clear that 
the poles of fl4"3| ) are approximately at c = ±c r — ij, where c r = and 7 <C c r ; both of the 
approximate poles are in the lower half plane, as expected, since we are considering a retarded 
propagator, which is required for the study of damping. 

To find expressions for c r and 7 we write H\ l \c) = c H\ l \c 2 ), so that, after the approxima- 
tion which follows, we do not lose the property that is an odd function of c. Now, since 
the functions and H\ are small compared to the regular term, we can make an approxi- 
mation similar to the one which leads to the well-known Breit-Wigner formula; we expand the 
quantities c 2 H^(c 2 ) and c 2 H\ l \c 2 ) which appear in (|43|) about c 2 = | keeping only the lowest 
order term. Thus we arrive at 



R 

S G (qo, |q|) ~ <( <l 2 - 



' 2(1 + J) 2#« ,2c H\ l) 



(44) 

3R 3R V3i? 1 

where the bar denotes quantities evaluated at c = 4=. ([|4]) exhibits the desired quadratic 
behaviour and has poles at c = ±c r — 27, where 



2(1 + J) 2HP _ H\ l) 

\ SR 3R 7 ~ y^3R 

and 7 2 has been neglected. c r and 7, as functions of temperature, are plotted in figfja and fig.|3|b 

respectively. They both have the expected behaviour, which reassures us about the validity of 

our approximation. In terms of c r and 7 we may recast ( fill ) as follows 

Sg(Qo, |q|) ~ I 7T2 ko - c r v f q 2 + 2i 7 Vf qo |q|l 1 (45) 
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or equivalently as 



Sa „ j 1 J— j (46) 

cv it |q| [ go - c r «f |q| + z 7 "^F |q| go + c r Vp |q| + «7 fF |q| J 

which manifestly has the correct form for a bosonic propagator with damping, the line-width 
having an explicit |q|-dependence. 
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Ill 



At this stage we have to point out that the alternative approximate form for Sg given 
1C] is incorrect, since it fails to exhibit the important property (traceable to ) that the 



imaginary part should be an odd function of q . The pole of Sq at c = — c r — 27 was overlooked 
in 



10f and instead a pole at c = — c r -M7 was used, which is actually on a Riemann sheet distant 
from the physical boundary This mistake had the effect of replacing the 2i 7 vp go |q| term in 
( |45| ) with something proportional to q 2 . Such a form then does, of course, have a simple local 



Fourier transform. In contrast ( |45l) does not but, as noted in the start of the section, the odd-go 
piece will integrate to zero in S e g. Damping can only be manifested in the equation of motion 
and in this case, when included, it will lead to a non-local term in the effective wave-operator 
(the transform of Sq 1 ). Nevertheless, we do have a simple and reliable expression ( pEpp for Sq 
in momentum space, which shows a particular form of damping and a predicted temperature- 
dependence for the line-shape. As noted in the Introduction, the non-local form of Sq 1 appears 
to present a fundamental difficulty for the inclusion of electromagnetism. 

We now turn to the treatment of the other terms on the right hand side of (^) . 

7.2 The remaining Landau terms involving and 

The fact that Sq is so sharply peaked around g ~ ±-^t>^|q| suggests that, as far as the 



dynamics of 6 described by (42) are concerned, it should be a good approximation to evaluate 
the c-even functions and at c 2 = |. The corresponding imaginary parts are very 
small, c-odd functions, as was, and we shall from now on simply ignore them. In this 



case all the terms in fl41|) become local and the Fourier transforms to x-space can all be done, 



leading to the effective Lagrangian 



Am 4 Am 



_jPp_ ( g 2 jy(l) ■ 2 g( W ) 2 ff(2) + [(V g ) 2 ] 2 ff(3)\ m 

+ 8mv F { r + Am r + (4m) 2 r )> 10 



where 



M\T) = ^ (1 + J) 



is flTTD as modified by the contribution from the Landau term, H)?> = H^\c 2 = |) (j = 1, 2, 3) 



and we have modified the 6 term as suggested after ([4lD . We note that the first part of (| 
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not involving the H)?> (j = 1,2,3) can be written as a TDNLSL just as in section 5, with a 
potential 



V 



(p - 2M 

w 



A2 



Choosing a convenient normalisation, the equation of motion for 9 is then 



0_ 

Of 

V. 



P 



2mvp 



(vey 

4m 



V0 

2m 



2mvp \ 



Am 



where 



p = 2M' - N[6 
It is tempting to interpret the quantity 



4m 



Pn 



2mvp 



4m 



(48) 



(49) 



(50) 



in (f|^) as an additional "density", not described by the NLSL. On this interpretation, we 
would expect the corresponding additional term inside the divergence in (|4"8"D to represent the 
associated "current" , which would require 

\2 



Pn~V r , 



2mv F \ 



4m 



2m 



(51) 



Equations flBTj) and (EI) are compatible, if 



V0 



2m* 

m* = (RP/HW)m 

h^hw = (spy. 



In this case, ([48]) becomes the continuity equation of the two-fluid model 



Jj(p + p n ) + V.(pv + p„v r , 







(52) 

(53) 
(54) 

(55) 



where p is the superfluid density and the superfluid velocity is v = V6*/2m as usual. 

The crucial requirement is (|54]). Fig.f|a shows the ratio Hf> / (H^ 2 ') 2 over the range of 
b we are considering, and it can be seen that (53) is quite well satisfied. Figfjb shows the ratio 
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I H^> as a function of b, from which we see that m* in fl53p varies quite rapidly with T at 
low T, and then approaches the constant value 5m. 

It is easy to show that there is a TDNLSL for the normal fluid from which the continuity 
equation for the normal fluid can be derived as an equation of motion. One can then use the 
"minimal coupling" prescription to take electromagnetism into account. 

In summary, we have argued that the effect of the Landau terms can be approximately 
incorporated by treating them as a constant, evaluating them at c 2 = 1/3. Considering the 
resulting effective Lagrangian, it seems to us to describe a two-fluid model in which the super- 
fluid component could be described by a TDNLSL with a temperature dependent potential, 
and the normal component by another TDNLSL with an appropriately chosen potential and a 
definite velocity predicted by fl52|). 
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Figure 1: (a): ^p- is very close to 1 for any temperature in the region which is considered here, 
b = y G [0,0.6]. (b): \J(T)\ is much smaller than 1 for any temperature in the same region. 
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Figure 2: (a): c 2 \H^\ is negligible compared to | ~ 0.333, for any temperature b £ [0, 0.6] and 
|c| < 1. For |c| > 1, it can be shown that the c-even function c 2 \H^\ is always much smaller 



than |. 
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Figure 2: (b): c 2 \H^\ is negligible compared to | ~ 0.333, for any temperature b G [0, 0.6] and 
\c\ < 1. For \c\ > 1, #f } = 0. 
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Figure 3: (a): c r is very close to the original location of the absolute value of the poles 
0.577, for any temperature in the considered region, (b): The damping 7 is negligible compared 
to the real part of the pole, for any temperature in the considered region. 
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Figure 4: (a): As the temperature rises approaches 1. (b): As the temperature 

rises H^/H^ = m* jm approaches 5. 
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